A new study of the polarized parton 
densities in the nucleon 



Elliot Leader 

Birkbeck College, University of London 
Malet Street, London WC1E 7HX, England 
E-mail: e.leader@physics. bbk. ac.uk 

Aleksander V. Sidorov 

Bogoliubov Theoretical Laboratory 
Joint Institute for Nuclear Research 
41980 Dubna, Russia 
E-mail: sidorov@thsunl.jinr.ru 

Dimiter B. Stamenov 

Institute for Nuclear Research and Nuclear Energy 
Bulgarian Academy of Sciences 
blvd. Tsarigradsko Chaussee 72, Sofia 1784, Bulgaria 
E-mail: ■ stamenov @inrne. bas. bg 

Abstract 

We present a new next-to-leading order QCD analysis of the world data on 
inclusive polarized deep inelastic lepton-nucleon scattering adding to the old set 
of data the final SMC results, the HERMES proton and very recent SLAC/E155 
deuteron data. We find an excellent fit to the data and present results for the 
polarized parton densities in different factorization schemes. These results are in 
a good agreement with what follows from the theory. We have also found that 
the main effect of the newly incorporated data is a better determination of the 
polarized gluon density. 



Recently JT|, [| we carried out a next-to-leading order (NLO) QCD analysis of the 
world data on polarized inclusive deep inelastic scattering (DIS). Since then the SMC 
group at CERN have re-analyzed their small-x data, using a totally new method of 
analysis, and published the final results of their experiment ||. At the same time new 
data have emerged from the HERMES Collaboration at DESY |4[] and very recently 
from the E155 experiment at SLAC f|. 

The SMC group have themselves carried out an NLO QCD analysis based on 
their final results and most of the world data available at that time (with the excep- 
tion of the final results from the SLAC/E143 experiment J7|). Given that there are 
significant changes in the SMC data, and the high precision of the new HERMES and 
SLAC/E155 data, we felt it necessary to redo our analysis. 

We have found that the newly incorporated data improve the determination of the 
polarized parton densities, in particular, of the poorly known gluon density. In this 
letter we present a brief summary of the general trends in our new results as well as a 
concise description of our present knowledge of the polarized parton densities. 

All details of our approach are given in |J. Here we simply recall that we carry 
out our analysis in three different factorization schemes:MS, AB (Adler-Bardeen) |§ 
and JET [|]. We then test the stability of our analysis by comparing densities in one 
scheme, as determined from the fit to the data, with these same densities evaluated 
using the transformation rules connecting densities in different schemes. 

We prefer to work with separate valence and sea parton densities and we check the 
stability of our results when we vary the parameter A corresponding to different flavour 
decompositions of the sea 

Au = Ad = XAs, (1) 

i.e., we confirm that the non-singlets Ag3, Ag§, the singlet AS and the gluon density 
AG are invariant. Note that this implies that the strange quark density 

As = i(AE-A 98 ) , (2) 
o 

is also invariant (does not change as A is varied) and it can be extracted from the data 
as well as AS. Of course, the valence quark densities Au v and Ad v are sensitive to the 
different assumptions about the sea. Nonetheless they are of interest for predicting the 
behaviour of other processes, e.g., polarized semi-inclusive DIS, polarized Drell-Yan 
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reactions, etc. 



For the input polarized parton densities at Ql = 1 GeV 2 we have adopted a very- 
simple parametrization 

Au v {x, Qq) = T]uA u X XUv^XjQq) , 

Ad v (x, Ql) = r] d A d x ad xd v (x,Ql) , 
ASea(x, Ql) = r]sAsx as xSea(x, Ql) , 
AG(x, Ql) = Vg A 9 x a *xG(x,Ql) , (3) 

where on RHS of (|3|) we have used the MRST unpolarized densities flQ |. The normal- 



ization factors Af in (§) are fixed such that rjf are the first moments of the polarized 
densities. 

The first moments of the valence quark densities rj u and rjd are fixed by the octet 
nucleon and hyperon f3 decay constants 



g A = F + D = 1.2573 ± 0.0028, a 8 = 3F — D = 0.579 ± 0.025 , (4) 
and in the case of SU(3) flavour symmetry of the sea (Am = Ad = As at Ql) 

Vu = 0.918 , Vd = -0.339 . (5) 
The rest of the parameters in (^) , 

{a u , a d , rjs, a s , % ,a g } , (6) 
have been determined from the best fit to the A±(x, Q 2 ) data. 

The numerical results of our fits to the world data on A±(x, Q 2 ) [3-5, 7, 12-15] are 
summarized in Table 1. The data used (161 experimental points) cover the following 
kinematic region: 

0.004 < x < 0.75, 1 < Q 2 < 72 GeV 2 . (7) 

As in our previous analyses P], Q the total (statistical and systematic) errors are 
taken into account. The results presented in Table 1 correspond to an SU(3) symmetric 
sea. Note that in this case a s = a,s and the first moment of the strange sea quarks, 
7fe = As(l,Qg) = 77 S /6. 

It is seen from the Table 1 that the values of x 2 /DOF coincide almost exactly 
in the different factorization schemes, which is a good indication of the stability of the 
analysis. The NLO QCD predictions are in a very good agreement with the presently 
available data on A± and g^, as is illustrated in the JET scheme fit in Figs. la-c. 
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Table 1. Results of the NLO QCD fits in the JET, AB and MS schemes to 
the world A± data (Qq = 1 GeV 2 ). The errors shown are total (statistical 
and systematic). 



Scheme 


JET 


AB 


MS 


DOF 


161 - 6 


161 - 


6 


161 - 6 


x 2 


128.3 


128.4 


129.4 


X 2 /DOF 


0.828 





.82* 


I 


0.835 




0.276 ± 0.030 


0.278 


± 


0.031 


0.257 ± 0.024 




0.077 ± 0.135 


0.065 


± 


0.136 


0.179 ± 0.107 


a§ 


1.524 ± 0.391 


1.668 


± 


0.403 


0.761 ± 0.188 


Vs 


- 0.032 ± 0.005 


- 0.029 


± 


0.006 


- 0.051 ± 0.006 


a g 


0.175 ± 0.452 


0.144 


± 


0.417 


3.6 ± 3.8 


% 


0.57 ± 0.14 


0.58 


± 


0.12 


0.07 ± 0.10 


AE(1) 


0.389 ± 0.037 


0.407 


± 


0.044 


0.275 ± 0.044 


a (l GeV 2 ) 


0.26 ± 0.05 


0.27 


± 


0.05 


0.28 ± 0.04 



Let us now comment briefly upon the state of our knowledge about the individual 
polarized parton densities. The extracted quark (valence and sea) and gluon densities 
at Qq — 1 GeV 2 are shown in Fig. 2 and Fig. 3, respectively. 

Au v (x, Q 2 ) and Ad v (x, Q 2 ) 

i) The valence quark densities are practically unchanged if the new (final SMC, 
HERMES proton and SLAC/E155 deuteron) data are incorporated in the analysis. 
(Note that Au v is very well determined once an assumption about the flavour decom- 
position of the sea is made.) 

ii) Au v and Ad v in the JET, AB and MS schemes coincide within the errors, in 
excellent agreement with what follows from the theory (they should be the same in the 
factorization schemes under consideration). 

iii) Although the first moment of Ad v (x, Qq), i]d, has been kept fixed in the analysis 
(see Eq. ([5])), the parameter a^, and therefore Ad v is not well determined from the 
present data (see Fig. 2). 

As(x,Q 2 ) and AE(x,Q 2 ) 

Note that these quantities are scheme dependent, 
i) As and AE are well determined from the data now (see Fig. 2). 
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ii) The first moments r] s in the JET and AB schemes are in a very good agreement. 
The same is valid for the first moments A£(l) (see Table 1). (We recall that according 
to the definition of the JET and AB schemes AE(1) as well as rjs should be the same in 
both schemes.) The corresponding densities As(x, Q 2 ) and AE(i, Q 2 ) in both schemes 
are slightly different because their higher moments are not equal. 

iii) The first moment of As in the MS scheme, (^s)ms' * s amios t twice as big as it 
is in the JET(AB) schemes, (?7s)jet(ab) : 

(Vs)m = -°- 051 ± °- 006 > (Vs)jet = -0.032 ± 0.005. (8) 

The polarized strange quark density is significantly different from zero indepen- 
dently of the factorization schemes used in the analysis. The result ([8]) demonstrates 
the extent to which the strange quark densities in different schemes can differ in the 
polarized case. 

Gluons AG(x, Q 2 ) 

i) The new set of data allows a better determination of the polarized gluon density 
in the JET and AB schemes. 

In previous studies \ 2 was quite insensitive to the value of a g in all factorization 
schemes under consideration, and the results of the fits corresponding to fixing a g = 0.6 
were presented in ]I|, £§. The results when Qjg IS £1 free parameter are summarized in 
Table 2. 

Table 2. The values of the parameters associated with the input 
polarized gluon density in the JET, AB and MS schemes for the 
new and old sets of data. 



Par am. 


Data set 


JET 


AB 


MS 


a g 


old 


1.8 ± 2.1 


1.7 ± 1.8 


1.6 ± 1.9 


Vg 


old 


0.24 ± 0.35 


0.28 ± 0.40 


0.30 ± 0.50 


a g 


new 


0.18 ± 0.45 


0.14 ± 0.42 


3.6 ± 3.8 


Vg 


new 


0.57 ± 0.14 


0.58 ± 0.12 


0.07 ± 0.10 



It is clear from the Table 2 that now both the value of a g and r? 9 , i.e., the shape 
and the normalization of the gluon density are much better constrained in the JET 
and AB schemes. The progress in extracting the polarized gluon density from the new 
data set is illustrated in Fig. 3 (JET scheme). 
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ii) In the MS scheme the gluons are still poorly constrained, x 2 continues to be 
insensitive to a large range of possible values for a g (see Table 2). The obtained value 
of rjg in this case is consistent with zero. However, if we fit the data with a fixed value 
of a g = 0.2 (in accord with the values obtained in the JET, AB schemes), we find 
rig = 0.62 + / — 0.42, in agreement with the values of this quantity in the JET and AB 
schemes, as expected theoretically. 

Axial charge a (Q 2 ) 

i) The good agreement between the values of a (Q 2 ) determined in the different 
schemes is demonstrated in Table 1, which illustrates how our analysis respects the 
scheme-independence of physical quantities. 

ii) The central values of a at Q 2 = 1 GeV 2 are slightly smaller than the correspond- 
ing values obtained from our previous analysis. 

We have also investigated the sensitivity of x 2 to the alternative possibilities of 
explanation of the spin of the nucleon via the spins of its constituents, namely: 
i) i] s ^ 0, f] g = , so that a AS << a 8 = 0.58 

ii) i] s = 0, f] g ^ , so that a « AS « a 8 = 0.58 

iii) ri- s ^ 0, r]g ^ 

The numerical results (JET scheme) are given in Table 3. 



Table 3. Sensitivity of x 2 to the absence of polarized strange 
quarks As or gluons AG at Q 2 = 1 GeV 2 . 



Assumption 


x 2 


- Vs 


Vg 


a 


Vs + 0, ri g = 


133.1 


0.065 ± 0.014 





0.19 ± 0.09 


Vs = 0, Vg ^ 


134.7 





1.33 ± 0.11 


0.27 ± 0.04 


Vs + 0, ^ ^ 


128.3 


0.032 ± 0.005 


0.57 ± 0.14 


0.26 ± 0.05 



It is seen from the Table that the present inclusive DIS data prefer small strange 
quark and small gluon polarized densities at Q 2 = 1 GeV 2 : 

jy- = -0.032 ± 0.005, % = 0.57 ± 0. 14. (9) 
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First moments of the spin structure functions 

Finally, we present our results for the first moments of the nucleon spin structure 
functions g± at Q 2 = 5 GeV 2 in the measured x range from 0.003 to 0.8 (AB scheme): 

0.133 for Proton 
0.039 for Deuteron (10) 
•0.048 for Neutron 

which are in excellent agreement with their experimental values || obtained from the 
world set of data (HERMES proton and E155 deuteron data not included): 

0.130 ± 0.003(stat) ± 0.005(syst) ± 0.004(ewo/) for Proton 
0.036 ± O.OOA(stat) ± 0.003(syst) ± 0.002(ewo/) for Deuteron (11) 
-0.054 ± 0.007(stat) ± 0.005{syst) ± 0.00A{evol) for Neutron 

The nucleon spin 

In the JET and AB schemes it is meaningful to interpret AE(1), the first moment 
of the singlet density AE(x, Q 2 ), as the contribution of the quark spins to the nucleon's 
spin. Our value of AS = 0.40 ± 0.04 is not far from the value 0.6 expected in certain 
quark models [JHJ. For the spin contribution from both quarks and gluons we have 
(JET scheme): 

-A£(l) + rj g = 0.76 ±0.14 (12) 

This value is consistent with 1/2 in two standard deviations. The more accurate 
determination of AS and rj g will answer the basic question how the spin of nucleon is 
divided up among the spin of quarks and gluons and their orbital angular momenta. 

In conclusion, we have re-analyzed the world data on inclusive polarized deep in- 
elastic lepton-nucleon scattering in NLO QCD adding to the old set of data the final 
SMC results, the new HERMES proton and SLAC/E155 deuteron data. As in the 
previous analysis it was demonstrated that the polarized DIS data are in an excel- 
lent agreement with the pQCD predictions for A±(x, Q 2 ) and gi(x,Q 2 ) in all the 
factorization schemes considered. The polarized parton densities have been extracted 
from the data. We have found that the main effect of the newly incorporated data is a 
better determination of the polarized gluon density, but in comparison with the other 
densities the uncertainty is still large. It follows from our analysis that the present 
inclusive DIS data prefer small strange quark and small gluon polarized densities, but 
they are significantly different from zero. 



/" dxg?(x,Q 2 = 5GeV 2 ) ftt = 

J0.003 
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To test more precisely the spin properties of QCD and to determine better the 
polarized densities, accurate data from both the inclusive (neutral and charged cur- 
rent) polarized DIS and the semi-inclusive processes in a larger kinematical region are 
needed. Finally, a direct measurement of the gluon polarization is necessary. We hope 
the COMPASS experiment at CERN, the future experiments at the hadron collider 
RHIC and the possibility of having a polarized proton beam at HERA will help to 
further a more profound study of the internal structure of the nucleon. 

This research was partly supported by a UK Royal Society Collaborative Grant, by 
the Russian Foundation for Basic Research, Grant No 99-01-00091 and by the Bulgar- 
ian science Foundation under Contract Ph 510. 
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Figure Captions 



Fig. 1. Comparison of our NLO results in the JET scheme for A^(x,Q 2 ) (a,b) and 
xgi(x,Q 2 ) (c) with the experimental data at the measured x and Q 2 values. Errors 
bars represent the total error. 

Fig. 2. Next-to-leading order polarized valence, strange and singlet quark distribu- 
tions at Q 2 = 1 GeV 2 with their error bands (JET scheme). The systematic errors 
are added quadratically. 

Fig. 3. Next-to-leading order polarized gluon density at Q 2 = 1 GeV 2 determined 
from the old (a) and the new (b) world sets of data (JET scheme). The error bands 
account for the statistical and systematic uncertainties. 
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Fig. la 




Fig. lb 




Fig. lc 



